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Abstract. We study quantum transport on finite discrete structures and we model 
the process by means of continuous-time quantum walks. A direct and effective 
comparison between quantum and classical walks can be attained based on the average 
displacement of the walker as a function of time. Indeed, a fast growth of the average 
displacement can be advantageously exploited to build up efficient search algorithms. 
By means of analytical and numerical investigations, we show that the finiteness and 
the inhomogeneity of the substrate jointly weaken the quantum walk performance. We 
further highlight the interplay between the quantum- walk dynamics and the underlying 
topology by studying the temporal evolution of the transfer probability distribution 
and the lower bound of long time averages. 
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1. Introduction 

Quantum walks (QWs) are attracting increasing attention in many research areas, 
ranging from solid-state physics to quantum computing [1]. In particular, QWs provide a 
model for quantum-mechanical transport processes on discrete structures; this includes, 
for instance, the coherent energy transfer of a qubit on an optical lattice [21 [3l lU [5] . The 
theoretical study of QWs is also encouraged by recent experimental implementations 
able to corroborate theoretical findings P, [3, E] . 

As in the classical random walk, quantum walks appear in a discrete [9] as in a 
continuous-time (CTQW) [lO] form; these forms, however, cannot be simply related 
to each other pjj. Now, standard CTQWs, on which we focus, can be obtained by 
identifying the Hamiltonian of the system with the classical transfer matrix which is, in 
turn, directly related to the Laplacian of the underlying structure. 

Another feature which CTQWs share with classical random walks consists in the 
strong interplay between the dynamics properties displayed by the walk and the topology 
of the substrate [121 [IS [H]- However, the dependencies turn out to be much more 
complex in the quantum-mechanical case: while the classical (simple) walk eventually 
loses memory of its starting site, the quantum walk exhibits, even in the asymptotic 
regime, transition probabilities which depend on the starting site. For this reason, often 
the parameters describing the transport are averaged over all initial sites, a procedure 
which allows a global characterization of the walk, while preserving its most important 
features. 

One of the quantities affected by topology is the mean square displacement of the 
walker up to time t. Classically, this quantity is monotonically increasing and depends 
(asymptotically) on time according to the power law (r^(t)) ~ t^. The value of the 
"diffusion exponent" f3 allows to distinguish between normal {(3 = 1) and anomalous 
{j3 7^ 1) diffusion [I5]. As for quantum transport, it is possible to introduce analogous 
exponents, characterizing the temporal spreading of a wave-packet [16]. However, even 
when they take place over the same structure, quantum and classical walks can exhibit 
dramatically different behaviours. In particular, the quantum wave propagation on 
regular, infinite lattices is ballistic, i.e. the root mean square displacement is linear 
in time. Such a quadratic speed-up of the mean square displacement is a well known 
phenomenon when dealing with tight-binding electron waves on periodic lattices [9] 
and, from a computational point of view, it constitutes an important feature since it 
could be advantageously exploited in quantum search algorithms [TTl [HI [19]. In the 
presence of disorder (either deterministic or stochastic) or finiteness, the sharp ballistic 
fronts are softened, a fact which may even lead to the localization of the quantum 
particle [201 [21] ■ It is therefore of both theoretical and practical interest to highlight how 
finiteness and inhomogeneity - often unavoidable in real systems - affect the particle's 
propagation. To this aim we analyze quantum transport on restricted geometries, where 
the restrictions arise from the (possibly joint) fractal dimension and finite extent of the 
substrate itself. By direct comparison with the classical case, we find that, on finite 
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substrates, the advantage of CTQWs is at short times only. Moreover, the lack of 
translational invariance weakens the CTQW performance, i.e. in such situations the 
average displacement increases more slowly with time. 

The finite discrete structures we consider and compare are the Dual Sierpinski 
Gasket (DSG), the Cayley Tree (CT) vide infra Sec. HI and the square lattice with 
periodic boundary conditions, i.e. the square torus (ST). These constitute representative 
topologies, providing examples of fractals with loops, of trees and of regular structures. 
The Dual Sierpinski Gasket will be treated in more detail; for this structure the 
eigenvalue spectrum of the Laplacian matrix is known exactly, allowing for some 
analytical estimates. Indeed, not only random walks, but also many dynamical 
properties of connected structures themselves (such as the vibrational structures and the 
relaxation modes) depend on the spectrum of their Laplacian matrix [22] . However, for 
CTQWs the set of eigenvectors also matters, which often makes analytical investigations 
cumbersome. 

It is worth underlining that focusing on discrete structures is not only suggested 
by solid-state applications: quantum computation is traditionally concerned with the 
manipulation of discrete systems. In particular, a discrete (and finite) state space makes 
the CTQW simulation by quantum computers, working with discrete registers, feasible 

mm- 

Our paper is structured as follows. After a brief summary of the main concepts 
and of the formula concerning CTQWs in Sec. [2l we describe the topology of the DSG 
in Sec. [3l Then, in Sec. HI we study the quantum-mechanical transport over the above 
mentioned structures, especially focusing on the average displacement and on the long 
time averages. Finally, in Sec. [5] we present our comments and conclusions. In the 
Appendix we derive analytical results concerning the average chemical displacements of 
CTQWs over hypercubic lattices, special cases being chains and square lattices. 



2. Continuous-time quantum walks on graphs 

Mathematically, a graph is specified by the pair {V, E} consisting of a nonempty, 
countable set of points V, joined pairwise by a set of links E. The cardinality of V 
provides the number A/" of sites making up the graph, i.e. its volume: \V\ = M . In the 
following, we focus mainly on finite graphs (A/" < oo) and we label each node with a 
lowercase letter i & V. 

From an algebraic point of view, a graph can be described by its adjacency matrix 
A, whose elements are 

A = h if (^'^■)^^ 

I otherwise. 

The connectivity of a node i can be calculated as a sum of matrix elements Zi = Aij. 
The Laplacian operator is then defined as L = Z — A, where Z is the diagonal matrix 
given by Zik = zAk- 



Dynamics of continuous-time quantum walks in restricted geometries 



4 



The Laplacian matrix L is symmetric and non-negative definite and it can therefore 
generate a probabihty conserving Markov process and define a unitary process as well. 
Otherwise stated, the Laplacian operator can work both as a classical transfer operator 
and as a tight-binding Hamiltonian of a quantum transport process [2H [25] . 

Indeed, the classical continuous-time random walk (CTRW) is described by the 
following Master equation [26] : 

d ^ 

■TiPkjit) =J2TkiPi,j{t), (1) 
1=1 

where pkj (t) is the conditional probability that the walker is on node k when it started 
from node j at time 0. If the walk is symmetric with a site- independent transmission 
rate 7, then the transfer matrix T is simply related to the Laplacian operator through 
T = -7L. 

Now the CTQW, the quantum-mechanical counterpart of the CTRW, is introduced 
by identifying the Hamiltonian of the system with the classical transfer matrix, H = — T 
[Tot [T3l [21] (in the following we will set h = 1). The set of states representing the 
walker localized at the node j, spans the whole accessible Hilbert space and also provides 
an orthonormal basis set. Therefore, the behaviour of the walker can be described by 
the transition amplitude ak,j{t) from state \j) to state \k), which obeys the following 
Schrodinger equation: 

d ^ 

-akj{t) = -tJ2Hkiai,,{t). (2) 

1=1 

If at the initial time t^ = only the state \j) is populated, then the formal solution to 
Eq. [2] can be written as 

akj{t) = {k\exp{-tm)\j), (3) 

whose squared magnitude provides the quantum mechanical transition probability 
7rfcj(t) = \akj{t)\'^. In general, it is convenient to introduce the orthonormal basis 
!'?/'„), n G [IjA/"] which diagonalizes T (and, clearly, also H); the correspondent set of 
eigenvalues is denoted by {\n}n=i,...,Ar- Thus, we can write 



(4) 



^{k\e-'^-'\^n){i^n\3) 

n=l 

Despite the apparent similarity between Eq. [1] and [2], some important differences 
are worth being recalled. 

First of all, the imaginary unit makes the time evolution operator U(t) = exp(— iHt) 
unitary, which prevents the quantum mechanical transition probability from having 
a definite limit as t 00. On the other hand, a particle performing a CTRW is 
asymptotically equally likely to be found on any site of the structure: the classical 
Pk,jit) admit a stationary distribution which is independent of initial and final sites, 
\imt^ooPk,j(t) = 1/Af. Hence, in order to compare classical long time probabilities 
with quantum mechanical ones, we rely on the long time average (LTA) p7] , defined in 
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Sec.lOl 

Moreover, the normalization conditions ioi pkj{t) and ak,j{t) read Yl^=iPk,j{t) = ^ind 

2.1. Average displacement 

The average displacement performed by a quantum walker until time t allows a 
straightforward comparison with the classical case; it is also more directly related to 
transport properties than the transfer probability 7ikj{t): It constitutes the expectation 
value of the distance reached by the particle after a time t and its time dependence 
provides information on how fast the particle propagates over the substrate. 

For CTQW (subscript q) starting at node j, we define the average (chemical) 
displacement {rj{t))q performed until time t as 

M 

k=l 

where i{k,j) is the chemical distance between the sites j and k, i.e. the length of the 
shortest path connecting j and k. We can average over all starting points to obtain 

(^it)), = JfJ2('Mr (6) 

i=i 

For fractals or hyperbranched structures it is more appropriate to use the chemical 
distance, rather than the Euclidean distance; for instance, the infinite CT (see Sec. H]) 
cannot be embedded in any lattice of finite dimension. For classical diffusion, it is well- 
known that the chemical and the Euclidean distances display analogous asymptotic laws 
for regular structures and for many deterministic fractals (e.g. the Sierpinski Gasket) 
^15j ; as discussed in the Appendix, this still holds for CTQWs on arbitrary (i-dimensional 
hypercubic lattices. 

For classical (subscript c) regular diffusion (on infinite lattices) the average 
displacement {r(t))c depends on time t according to 

{r{t))c - t'/'. (7) 

More generally, for scaling (fractal) structures we can define the so-called chemical 
diffusion exponent df^ and get ^15j: 

(r(t)),~ti/'^^. (8) 

Finite systems require corrections to these laws: for them {r{t))c does not grow 
indefinitely, but it saturates to a maximum value Vc |28j . 

2.2. Return Probability 

As it is well known, for a diffusive particle the probability to return to the starting 
point is topology sensitive, and it can indeed be used to extract information about the 
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underlying structure It is therefore interesting to compare the classical return 

probabihty Pk,k{t) with the quantum mechanical vrfc fc(t) (see also [291 [30]). One has 

Pk,kit) = (A;|exp(Tt)|A;) = ^ |(A;|^„)|'exp(-7tA„) (9) 



n=l 



and 



7rfc,fc(t) = |afc,fc(i)r 



\ {k\ipn)f exp{-i-ftXr. 



n=l 



(10) 



In order to get a global information about the likelihood to be (return or stay) at the 
origin, independent of the starting site, we average over all sites of the graph, obtaining 

Af , Af 



and 



Pit) 



nit) 



k=l 



n=l 



1 ^ 

rj=^'^k,kii) 



Af 



k=l 



1 

Af 



Af 
n,m=l 



Af 



■i7(-^n— Am)t 



J2m^n)\'m^m)f. (12) 



k=l 



For finite substrates, the classical pit) decays monotonically to the equipartition 
limit, and it only depends on the eigenvalues of T. On the other hand, 7f(t) depends 
explicitly on the eigenvectors of H [291130]. By means of the Cauchy-Schwarz inequality 
we can obtain a lower bound for fcit) which does not depend on the eigenvectors [30[[3T]: 



nit)> 



1 ^ 

r?^oik,kit) 



Af 



k=l 



1 

AP 



Af 
m.n=l 



-i-yiX„-Xm)t 



(13) 



Notice that Eqs. [TT] and [12] can serve as measures of the efficiency of the transport 
process performed by CTRW and CTQW, respectively. In fact, the faster pit) decreases 
towards its asymptotic value, the more efficient the transport. Analogously, a more 
rapid decay of the envelope of 7f(t) (or of |Q;(t)p) implies a faster delocalization of the 
quantum walker over the graph. By the way, we recall that, for a large variety of graphs 
[30] , the classical average return probability scales as pit) ~ while the envelope 
of |a(t)p, namely env [|a()f:)p], scales like /i being a proper parameter related for 
fractals to the spectral density. 

As can be inferred by comparing Eqs. [11] and [121, for quantum transport processes 
the degeneracy of the eigenvalues plays an important role, as the differences between 
eigenvalues determine the temporal behaviour, while for classical transport the long 
time behaviour is dominated by the smallest eigenvalue. Situations in which only a few, 
highly degenerate eigenvalues are present are related to slow CTQW dynamics, while 
when all eigenvalues are non-degenerate the transport turns out to be efficient [291 ]30] . 




Figure 1. Dual transformation from Sierpinski Gasket to Dual Sierpinski Gasket of 
generation g = 3. 



3. Dual Sierpinski gasket: Topology and eigenvalue spectrum 

Before turning to the dynamics of CTQW (and CTRW) on exemplary structures, which 
allow us to highlight the importance of inhomogeneities, some remarks on the spectra 
of the DSG are in order. 

The dual Sierpinski gasket is an exactly-decimable fractal which is directly related, 
through a dual transformation, to the Sierpinski gasket (SG). The DSG of generation g 
can be constructed by replacing each small triangle belonging to the SG with a node and 
by connecting such nodes whenever the relevant triangles share a vertex in the original 
gasket (see Fig. [1]). It is straightforward to verify that the number of nodes at any given 
generation g is N" = 3^. 

The dual transformation does not conserve the coordination number of the inner 
nodes (which decreases from 4 to 3), but it does conserve the fractal dimension df and 
the spectral dimension d, which are therefore the same as for the original Sierpinski 
gasket df = In 3/ In 2 = 1.58496... and J = 2 In 3/ In 5 = 1.36521.... 

As mentioned above, the knowledge of the eigenvalue spectrum is sufficient for the 
calculation of several interesting quantities concerning the dynamics of CTQWs. In 
general, any (finite) Hamiltonian H can be (at least numerically) diagonalized in order 
to obtain its spectrum. However, as the size of H gets large, the procedure gets to be 
time consuming and the precise numerical diagonalization may not be easy to perform. 
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Remarkably, the eigenvalue spectrum of the DSG Laplacian matrix can be determined 
at any generation through the following iterative procedure; for more details we refer to 
[311 EH] : At any given generation g the spectrum includes the non-degenerate eigenvalue 
A_/v' = 0, the eigenvalue 3 with degeneracy (3^~^ + 3)/2 and the eigenvalue 5 with 
degeneracy (3^~^ — l)/2. Moreover, given the eigenvalue spectrum at generation g — 1, 
each non- vanishing eigenvalue Ag_i corresponds to two new eigenvalues according to 

^ ^"^^P^ ; (14) 

both A^ and A^ inherit the degeneracy of Xg-i- The eigenvalue spectra is therefore 
bounded in [0,5]. As explained in [51], at any generation g, we can calculate the 
degeneracy of each distinct eigenvalue: apart from A^v" whose degeneracy is 1, there are 
2^ distinct eigenvalues, each with degeneracy (3^^''^^ + 3)/2, being r = 0, 1, ...,g — l, and 
2^ distinct eigenvalues, each with degeneracy (3^^^'^^ — l)/2, being r = 0, 1, ...,g — 2. 
As can be easily verified, the degeneracies sum up to A/" = 3^. Finally, notice that the 
distribution of eigenvalues and their degeneracies are non-uniform and that the spectrum 
is mult if r act al |34j . 



4. CTQWs on restricted geometries 

4-1. Transfer probability 

Results for the exact transition probability distribution vr^ for STs and CTs have 
already been given in Refs. [131 [311 132], where it was shown that i^kjif) depends 
significantly on the starting node. Results for ultrametric structures are given in |33j . 

It is worth recalling here that the Cayley tree (CT) can be built by starting from 
one node (root) connected to z nodes, which constitute the first shell. Each node of the 
first shell is then connected to 2; — 1 new nodes, which constitute the second shell and 
so forth, iteratively. Therefore, the M-th shell contains z{z — 1)^-^ nodes which are 
at a chemical distance M from the root. Thus, the CT is a z-regular loop-free graph. 
The numbers of sites in a CT of M shells is Mm = [z{z - 1)^ - 2]/{z - 2), hence the 
correlated fractal dimension log(A/M)/ log(M) goes to infinity for M — > 00, precluding 
the possibility of embedding very large CT in any previously specified Euclidean lattice. 
In the following we focus on finite 3-Cayley trees, which means that z is fixed and equal 
to three for any internal site of the graph; furthermore, the number of shells (also called 
generation) is finite (and therefore also the number of nodes is itself finite). 

In Fig. [2] we show our results for a DSG of generation g = 3 and we focus on 
the set of pairs given by (fn,l), where f„ denotes any of the two corners of the 
gasket of the n-th generation, with n ^ g (i.e., according to the labeling of Fig. [21 
Vq = l,Vi = 2, t>2 = 5,f3 = 14). Now, due to the symmetry the DSG is endowed with, 
for CTQWs starting from a given vertex, say the apex, the left and right corners are 
equivalent. As expected, iTkjit) does not converge to any definite value, but it displays 
oscillations whose amplitudes and average values get smaller as the distance between 
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Figure 2. Exact probability TTk.i{t) for the CTQW starting from the apex (site j = 1) 
to reach sites fc = 1, 2, 5 and 14. The fc-sites are the left corners at generations g — 0,1,2 
and 3, respectively. 

the sites 1 and f„ increases. This suggests, at least when starting from a main vertex, 
that the CTQW stays mainly localized at the origin and its neighbourhood. 

We corroborate this by looking at the temporal evolution of tt^ for the DSG and 
by comparing it to the TTk,j{t) pertaining to the CT and the ST of comparable size A/". 
Figures [3l H] and [5] show 3D pictures of Hkj (t) at different moments t (belonging to the 
short-time regime). On both the x and the y-axis, k and j label the nodes of the graph 
in such a way that at the point {k,j) on the xy plane the value of 7ik,j{'t) is presented. 
At the initial time t = 0, the transition probability 'n'k,j{t) is non- vanishing only on the 
diagonal, i.e. one has 7rfcj(0) = 6k,j; at later times, 'n'k,j(t) spreads out non-uniformly, 
according to the topology of the substrate. In particular, for the DSG (Fig. [3]), a large 
fraction of vrfcj(t) stays in a region connected by bonds to the initial nodes; several 
peaks can be distinguished, whose heights decrease as the chemical distance between 
the pertaining sites gets larger. For the CT, the pattern representing TTk,j(t) is even 
more inhomogenous; as can be inferred from Fig. HI a quantum particle on the CT is 
located with very high probability on its initial node (except when starting from the 
central node) and, for the time scale considered, it is very unlikely to reach nodes outside 
its starting branch. On the other hand, for the ST (Fig. [5]) we notice that the spread 
of TTk,j{t) is rapid and regular: apart from possibly partial revival phenomena (see for 
example the snapshot for t = 87"^), the pattern for 7!'k,j(t) exhibits very low peaks. 
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Figure 3. Snapshots of the transition probabihty tt^ j (t) at different times t — 0,1,3, A 
(in units of 7^^) for the Dual Sierpinski Gasket of generation g = 3 {Af = 27). 



Indeed, peaks in 7f k,j{t) are a consequence of the constructive interference stemming 
from reflections at peripheral sites or (in the case of the torus) from the superposition 
of traveling waves which have crossed the whole (finite) graph. 

Finally, Figs. [3l H] and [5] also highlight the symmetry characterizing the quantum 
transfer probability, namely that TTj,k{t) = TTk,j{t), at all times. This can be derived 
directly from Eq. [3l recalling that H is itself symmetric and real. An analogous 
symmetry also characterizes the classical distribution Pj,k{t) for all the cases analyzed 
here. 

4-2. Average displacement 

The dynamics of quantum particles on non-regular structures has been investigated 
in several works meant to analyze the quantum dynamics of tight-biding electrons 
in quasicrystals, in aperiodic and quasi-periodic chains and in random environments 
[20I [2TI [36] . There, the highlighted dramatic deviations from the ballistic behaviour 
(expected for regular, infinite lattices) range from anomalous to super diffusion, to 
decoherence, and even to Anderson localization [37] . 



Dynamics of continuous-time quantum walks in restricted geometries 



11 




Figure 4. Snapshots of the transition probabihty iTk.j (t) at different times t — 0,1,3, A 
for a Cayley Tree of generation g ^ 3 {Af = 3 x 2^ — 2 = 22); time is give in units of 
7~^. Notice that the distribution is locahzed on special couples of nearest-neighbours 
sites and that also reflection effects appears. 

Here, we consider the case in which non-regularity stems from the intrinsic spatial 
inhomogeneity of the substrate (for DSGs and CTs) and we also study quantum walks 
on STs which allow to evidence the role of finiteness. From the experimental side, the 
importance of such factors (spatial inhomogeneities and finiteness of the sample) has 
been increasingly recognized (see e.g. p8]), so that it is of great interest to understand 
to what extent quantum transport is influenced by them. 

First of all, we notice that the fact that vrfc j(t) does not attain a stationary 
distribution also causes the average displacement {rj(t))q not to necessarily increase 
monotonically with t. Moreover, due to reflection effects, we expect the mean value 
(r(t))g to overestimate the displacement performed by a CTQW that started from a 
peripheral site. Indeed, one finds for the DSG and the CT that {r{t))^ is larger than 
{rj{t))q, j being any corner of the gasket (see Fig. [6]) or any peripheral site, respectively. 
As for the ST, {r(t))^ trivially equals {rj(t))q, for all j. 

As mentioned above, for classical diffusion the average displacement grows 
continuously from zero to a maximum value fc which, due to equipartition, is just 
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Figure 5. Snapshots of the transition probabihty iTk.j (t) at different times t ~ 0,1,3,4 
(in units of 7^^) for the Square Torus of Unear size i = 5 {JV ~ 25). Notice that the 
distribution spreads very rapidly and regularly over the whole structure. 



the mean distance among sites: 



1 ^ 



J] 



(15) 



k,j=i 



Despite of the oscillating behaviour of {r{t))^, we can obtain an analogous constant 
value Vg, around which the average displacement eventually fluctuates, which reads 

hm ^ I dt^ (16) 



Otherwise stated, {r{t))^ eventually reaches a "stationary regime" in which it fluctuates 
around a constant value (see Fig. [7]). 

Of course, and Tc depend on both the topology and the size of the substrate, 
and they diverge as TV ^ 00. From the remarks of Sec. 14. H we expect that quantum 
interference arising from reflection affects Vg, making it smaller than Tc. Indeed, for 
the CT and the DSG, Fig. [7] clearly shows that < r^, this is especially apparent 
for the CT where ~ 8.9 (calculated from Eq. [T5|) is nearly four times larger than 
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'0 2 4 6 8 10 12 

Figure 6. Average displacement {rj{t))q for a quantum walker which started at the 
j-th site on a DSG of generation <? = 5. The main figure focuses on short times, while 
in the inset a wider temporal range is considered. Different colours and thicknesses 
distinguish different initial sites j, as shown in the legend; the labeling is the same as 
in Fig. [21 Notice the appearance of local minima from t « 7^^ onwards. 

rg 2.4. Conversely, for the ST, where interference only stems from the superposition 
of waves which have crossed the whole substrate, we find that Vc and Vg are eventually 
comparable. Therefore, we expect that on structures endowed with reflecting boundaries 
(i.e. peripheral nodes of low connectivity), at sufficiently long times, the expectation 
value of the average distance reached by a quantum particle is strictly smaller than the 
average distance Vc among the nodes. 

The classical and the quantum cases are further compared in Fig.[8]which shows the 
ratio {f'{t))^/{r{t))^: one can notice that, for significant times (t > l7~^), the classical 
average displacement is strictly lower than the quantum-mechanical one up to time 
t ^ 47-^^ ^ 97-1 and t ^ 187-^ for CT, ST and DSG, respectively. 

Thus, we can conclude that, on restricted geometries such as those analyzed here, 
CTQWs can spread faster than their classical counterpart, although the advantage is 
significant only at relatively short times. Moreover, the spatial homogeneity enhances 
the speed-up; especially for CTs and, in general, for tree-like structures, the large number 
of peripheral sites gives rise to localization effects which significantly reduce Vg. 

Finally, we stress that analytical results on the average displacement performed 
by a quantum particle on discrete structures are rather sparse (see e.g. [161 EH]); in 
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5 10 15 20 

Figure 7. Classical (dashed line) and quantum (continuous line) average displacement 
for the DSG (.9 = 5, TV = 243), the CT (5 = 6, = 190) and the ST (L = 16, 
J\f = 256). The main figure focuses on the short-time regime, while the inset also 
shows the long-time regime. 

the Appendix we prove that on infinite c?- dimensional hypercubic lattices both the 
average chemical displacement defined in Sec. 12.11 and the Euclidean displacement 
depend linearly on time and that this kind of behaviour survives, at short times, also 
for finite lattices. 

4-3. Average return probability 

The average displacement for CTQWs already highlighted some aspects of the role of 
inhomogeneities for transport processes. Now, we obtain further insights by considering 
the average return probability. 

For the DSG we can get p{t) without numerically diagonalizing L, since it only 
depends on eigenvalues which can be calculated iteratively. Figure [9] displays the 
averaged probabilities p{t), Tf{t) and |a;(t)p - numerically evaluated from Eqs. [HI [12] 
and [131 respectively - as a function of time, obtained ioi g = 5. The classical p{t) 
decays monotonically to the equipartition value l/A/", while the quantum- mechanical 
probabilities eventually oscillate around the value 0.7, which is larger than 3~^. 
Although the amplitude of fiuctuations exhibited by the lower bound is larger than 
that of the exact value, the agreement between the two quantities is very good. In 




Figure 8. Ratio between classical and quantum average displacement for the DSG 
{g — 5), the CT [g = 6) and the ST (L = 16), as shown by the legend. The main 
figure focuses on the short-time, while the inset also shows the long-time regime. In 
particular, for the square torus, the ratio {r{t)) ^/ {r{t)) ^ is at first smaller than unity 
and then it oscillates around unity; the highest peaks arc signs of (partial) revivals. 



particular, the positions of the extrema practically coincide and the maxima of vf (t) are 
well reproduced by the lower bound. An analogous behaviour was found also for other 
graphs, such as square lattices ^13j, Cayley trees [21] and stars [29j. Notice, however, 
that for the square lattices the lower bound turns out to be exact while for Cayley trees 
and for stars it is only an approximation, which, moreover, turns out to be less accurate 
than what we find here for the DSG. 

On short times {t < 87^^) it is possible to construct the envelope of Tf{t), which 
depends algebraically on t. The exponent is ~ —0.82, to be possibly compared with 
d/2 ^ —0.68 which is the exponent expected classically for the infinite DSG. The decay 
of the average return probability 7t{t) for the ST can be estimated as well: its envelope 
goes like (classically as p(t) ~ t~^)[l3l [30], implying a faster delocalization of the 
CTQW over the graph. 

Interestingly, for the DSG, the overall shape of if{t) does not depend significantly 
on the size of the gasket (Fig. [TU]) . In fact, the behaviour of 7t{t) is mainly controlled by 
the most highly degenerate eigenvalues. These do not change when increasing the fractal 
size (i.e. its generation). These values are: 3 with degeneracy mg{3) = (3^^^ -|- 3)/2, 5 
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Figure 9. Average return probability 7r{t) for the DSG of generation g = 5 on a log- 
fog scale. The comparison with the classical p{t) evidences that the classical random 
walk spreads more efficiently than its quantum- mechanical counterpart. The dashed 
line represents the envelope of 7f(i). 



with degeneracy mg{5) = (3^ ^ — l)/2, and (5 ± ^/13)/2 with degeneracy mg_i(3), see 
Sec. [31 

4.4- Long time averages 

As underhned in Sec. |2l the unitary time evolution does not allow a definite long-time 
limit for vr^ Then, in order to obtain information about the overall spreading of 
quantum walks, it is advantageous to use the long time average (LTA): 

XA:,i= lim- / dtTlkjit) =y2Sx,„X,„{k\iJn){i^n\j){j\i^m){'ipm\k), (17) 
" 71, m 

where 5a„,a™ equals 1 for A„ = and is zero otherwise. The LTA of vf (t) follows as 

^^T^oof r^^^W = K,xJ{Mn)\' (18) 

^ n,m,j 

for which we obtain a lower bound which does not depend on the eigenvectors 



1 . 



n,m 
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Figure 10. Average return probability 7f(t) for tlie DSG of generation g = 4 (bright 
colour) and g = 5 (dark colour). Its lower bound |a(i)p (dashed line) and the classical 
p{t) (dotted line) are also depicted, as shown by the legend. 



We first consider the DSG for which Fig. [TT] shows Xkj as a contour plot, whose 
axes are labeled by the nodes k = 1, ...,Af and j = I, Bright colours correspond 

to large, dark ones to small LTAs. First of all, we notice that the LTAs are far from 
being homogeneous and, hence, are not equipartioned. In particular, the values on the 
main diagonal are high, meaning that CTQWs have a high LTA probability to be at 
the starting node. 

The inhomogeneity of the pattern mirrors the lack of translation invariance of the 
DSG itself. For instance, v being the label assigned to any vertex of the main triangle, 
Xv,v is a global maximum; off-diagonal local maxima correspond to couples of connected 
nodes belonging to different minor triangles of generation g — 1. This allows to establish 
a mapping between the pattern of Xkj and the structure of the relevant DSG. Indeed, as 
suggested by the white delimiting lines in Fig. [11], the patterns of the LTA distributions 
exhibit self-similarity. 

As for X and its lower bound xib, we recall that the former can be calculated 
numerically, once all eigenvalues and eigenvectors of the Laplacian operator are known 
(Eq. [T8l) . while for the latter the knowledge of the eigenvalue spectrum is sufficient 
(Eq. [19]). Since the spectrum of the DSG is known, we can calculate xw analytically. 
Recalhng the results of Sec. [3] at generation g the spectrum of L displays Af distinct 
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Figure 11. Limiting probabilities for the DSG of generation 5 = 3 (top) and g — 4: 
(bottom), whose volumes are Af — 27 and Af — 81, respectively. The white lines enclose 
the limiting distributions for gaskets of smaller generations. Notice that the global 
maxima lay on the main diagonal and correspond to j = 1, 14, 27 and to j = 1, 41, 81 
for g ~ 3 and for (7 = 4, respectively. 
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eigenvalues, where 



9-1 9-2 

A)' = ^2" + ^2" + l = 3x 2'^-^ - 1. 

r=0 r=0 



We call the set of distinct eigenvalues {Xi}^^^ j^. Being m(Aj) the degeneracy of the 
eigenvalue Aj, we can write 



n,m=l 



Af 



n=l 



i=l 



Now, we go over to the space of distinct degeneracies, each corresponding to a number 
p of distinct eigenvalues and we get the final, explicit formula 

1 

X>Xib = jn ^[m(r)]V("^(r)) 



r=0 





[E 

I r=0 






2 


1 

3^ 





+ 3' 



ff-2 
r=0 



-39-r-l _ 1" 




X 2'' + ij 


2 





> 



Interestingly, in the limit g 



7 2_ 
00, the LTA X is finite: 



39 



(20) 



X > lim xib 



1 

14' 



and Xife reaches this asymptotic value from above. 

In Fig. [12] we show, as functions of A/", x ^"^^ its lower bound, calculated from 
Eq. [TS]and from Eq. [201 For comparison, the same quantities obtained for CTs and STs 
are also depicted. In the latter case, due to the regularity and periodicity of the lattice 
[ID] , the lower bound actually coincides with the exact value. For all cases considered, 
X is larger than the equipartion value (given by the dashed line). 

The inset of Fig. [T2] shows the ratio 



X_ 
Xib 



Obviously, the closer rj is to 1, the better xit approximates x- this sense, the lower 
bound calculated for CTs is not as good an approximation to x ^s it is for the DSG 
and for the ST. For the CT, xib definitely underestimates, being about half the exact 
value of X- The quantity r]{Af) may act as a measure of the inhomogeneity of a given 
substrate. Practically, when dealing with a large sized, sufficiently regular structure, we 
can get information about the localization of a quantum particle moving on it simply 
through xib, thus avoiding the (lengthy) evaluation of the eigenvector set. 
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Figure 12. Long time average x and its lower bound xib calculated according to 
Eq. [20] for the Dual Sierpinski Gasket (black diamonds and stars) , the Cayley Tree 
(red triangles and crosses) and the Square Torus (blue crosses and open circles) versus 
the size of the structure Af; lines are guides for the eye. The dashed line represents 
the equipartition value 1/Af. Inset: ratio 77 as a function of JV. Notice that r]{AJ') = 1 
holds not only for the periodic square lattice, but for all hypercubic lattices. 



5. Conclusions 



We investigated the behaviour of continuous-time quantum walks (CTQWs) on finite 
discrete structures characterized by different topologies; we considered the Square Torus 
(ST), the Cayley Tree (CT) and the Dual Sierpinski Gasket (DSG). 

The interplay between the quantum-walk dynamics and the underlying topology 
was deepened by studying, in particular, the temporal evolution of the transfer 
probability distribution and the ratio x/Xib as a function of the substrate size. The 
latter turns out to be significantly sensitive to the inhomogeneity of the substrate, 
from which we can infer that lower bound estimates are especially reliable for regular 
structures. 

From an applied, as well as theoretical, perspective, the average displacement of 
the walker, as a function of time, also plays an important role. This quantity is not only 
directly related to the transport properties, but it also provides information about how 
fast the walk explores the underlying structure, allowing an immediate comparison with 
the classical case. We found that at short times, CTQWs can spread faster than their 
classical counterparts, although spatial inhomogeneities and finiteness jointly reduce this 
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effect. In the Appendix we prove that for infinite (i- dimensional hypercubic lattices, at 
long times both the average chemical and the Euclidean displacements depend linearly 
on time (i.e. the motion is ballistic); for finite lattices this kind of behaviour holds at 
relatively short times only. 
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Appendix A. Average chemical displacement on hypercubic lattices 

Here we consider infinite (i-dimensional hypercubic lattices and, by exploiting their 
translational invariance, we prove that on them the average chemical displacements of 
CTQWs, as defined in Sec. 12. H depend linearly on time. We first focus on the infinite 
discrete chain, then we consider the generic c?- dimensional case and finally we analyze 
the two-dimensional lattice. 

For a ring of length A/", by exploiting the Bloch states, we have [H] 



/AT ^ 

where A; is the l-th eigenvalue of the Laplacian matrix L associated with the ring. In 
the limit J\f —>■ oo we are allowed to replace the sum over / by an integral, obtaining 



lim akj{t) = ^'=-^e-*^V,_,(2t), (A.2) 



where Jk{z) is the Bessel function of the first kind. In the calculation of the transfer 
probability T^kjif) the phase factor vanishes and we have Tikjit) = Jk-jC^t): which can 
be restated as 

7r,.,o(t) = 4'(2t), (A.3) 

due to the translational invariance of the structure. Clearly (in agreement with TTkfiit) 
being a probability distribution), one has for all t 

+00 +00 oo 

E vr,,o(t) = Yl ^'(2t) = Jo'(2t) + 2^ J2(2t) = 1, (A.4) 

fe=— oo A:=— oo k=l 

\k 



the last equality being based on Jk{z) = {—l) J^k{z) and on Eq. 8.536.3 in [12]. 

Now, the average chemical displacement of a CTQW which starts from and moves 
on an infinite chain (subscript q, 1) follows from Eq. [5] as 

(ro(t)),,i = (r(t)),,i = ^ i{k, 0) Jl{2t) (A.5) 

keV 



J2 \k\Jl{2t) = 2Y,kJl{2t). 
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Here, in the first equality we dropped the subscript due to the equivalence between 
the sites and in the last equality we exploited the symmetry of the Bessel functions, 
J-ki^) = Jki^)- Now, recalling the recursion formula Eq. 8.471.1 in [12] 



2k 

Jk-i{z) + Jk+i{z) = — Jkiz), (A. 6) 

we can write 

oo oo 

2j2kJl{z)=zJ2[Jk-i{z)Jk{z) + Jk{z)Jk+i{z)] = zJ{z), (A.7) 

k=l k=l 

by defining the function J^{z). Hence 

{rk{^)),,i=zJiz), (A.8) 

where we put 2t = z. The squared Bessel function J'^{z) is almost everywhere positive 
and the analysis of its zeros allows to state that, for t > 0, the sum appearing in the 
left-hand-side of Eq. IXISl is strictly positive; the same holds therefore for J'{z), for which 
we also notice from Eq. IA.7I that J^{0) = 0. Moreover, through the following recursion 
formula, Eq. 8.471.2 in 



2^Jfc(2;) = Jk-i{z) - Jk+i{z), (A.9) 
it follows by directly differentiating J^{z) and rearranging the terms 

dz 2 z ' 

where in the last expression we again used Eq. IA.6l for k = 1. The indefinite integral of 
Eq.Oiis (see Eq. 5.53 in [42J) 

J{z) = z Jl{z) + z Jl{z) - Uz).h{z)+C, (A.ll) 

as can be simply verified by differentiating Eq. lA.llI and using Eqs. IA.6I and IA.9I 
Furthermore, since i7(0) = 0, we have C = 0. 

Therefore, the following, for us fundamental, relation holds: 



Y,kJl{z) = ^-[zj',iz) + zj',iz)-Mz)Mz)]. (A.12) 

k=l 

Now, from Eqs. IA.5l and lA.12l we get the exact expression for the average chemical 
displacement 

(r = z [z Mzf + z J^{zf - Mz)Ji{z)\ . (A.13) 

For large z = 2t (i.e. long times) we can use the expansion (see Eq. 8.451.1 in [12] ) 



2 

Jk(z] 

TXZ 




(A.14) 



Consequently, inserting Eq. IA.14I for Jo{z) and Ji{z) into IA.131 we infer that the long 
time behaviour of the average chemical CTQW displacement on an infinite chain obeys 

4f 

(r(t)),,i - -. (A.15) 

71 
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Figure Al. Average chemical displacement {r{t))^ for the torus calculated according 
to Eq. [6] (symbols); the line represents the Euclidean average displacement for L = 19. 
The dotted and the dashed lines highlight the linear dependence on t exhibited by the 
average chemical distance and by the average Euclidean distance, respectively. 



This result is consistent with findings reported in |39j for the average square 
displacement. 

Let us now consider higher dimensional hypercubic lattices. Again, without loss of 
generality, we can assume the CTQW to start from the point = (0,0,...,0) so that 
the chemical distance attained by a walker being at the generic site k = {ki,k2, kd) 
is £(k, 0) = \ki\ + \k2\ + ... + \kd\. Furthermore, on a hypercubic lattice, assuming 
symmetric conditions in all directions, the probability distribution 7rko(i^) factorizes 
into the d-independent one-dimensional distributions Hk^flit): 

d 

vrk,oW = T^k^fi{t)T^k2fl{i) ■■■T^ka,o{t) = Yl'^k.fii't)- (A.16) 

i=i 

Hence 

d d d 

(ro(t)),,. = iY. \k,\),,d = = = rf(ro(t)),,i. (A.17) 

j=i j=i j=i 

In the last relation we used the fact that = since for each j only the 

distribution vr^^ o(^) matters, the other distributions adding up to a factor of unity each. 
Hence 

{r{t)),,d = d{r{t)),,,r^^. (A.18) 
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In particular, for the square lattice we have 

(r(t)),,2 ~ -, (A.19) 

71 

which was used in Fig. lAll (dashed line) to fit data relevant to the average chemical 
displacement performed by a CTQW on square tori of different (finite) sizes. As can 
be seen from the figure, the ballistic behaviour also holds for finite lattices, but for 
relatively short times only: at longer times the finiteness of the lattice starts to matter 
and the product of Bessel functions in Eq. IA.16I ceases to be a good approximation of 
the transfer probability. When the waves associated with CTQWs have crossed the 
whole lattice, interference effects start to occur and {r{t))q^2 exhibits a non-monotonic 
behaviour. From the same figure we also notice that the 0{l/t) contributions of Eq. IA.14I 
get to be negligible for t > 1 7"^. 

In Fig. lAll we also show data for the average Euclidean displacement which displays 
a ballistic behaviour at short times as well. Indeed, for a hypercubic lattice of arbitrary 
dimension d, the following relation holds (see e.g. [l3j) 

i=^(k,j)<||k-j||<£(k,j) (A.20) 

where ||k — j|| denotes the Euclidean distance between the lattice points k and j chosen 
arbitrarily. By averaging each term of the previous equation with respect to the transfer 
probability 7rkj(t) (we can again exploit the translational invariance of the substrate and 
fix j = 0), we find that the average Euclidean distance also scales linearly with time with 
a multiplicative factor bounded between 4:d/{\/3Ti) ~ 2.31^/7? and Ad/n. In particular, 
for the square torus of size L = 19 considered in Fig. lAll we find that at relatively short 
times the average Euclidean distance scales as Gt/ir. 
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